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Abstract 

We review several approaches to define and quantify decoherence. We find 
that a measure based on a norm of deviation of the density matrix is appro- 
priate for quantifying decoherence for quantum registers. For a semiconductor 
double quantum dot charge qubit, evaluation of this measure is presented. For 
a general class of decoherence processes, including those occurring in semicon- 
ductor qubits, we establish that this measure is additive: It scales linearly with 
the number of qubits in the quantum register. 



1 Introduction 

Decoherence PQHHS1I31IS1IS1IZII51IS1 0D1 IHl H21 US Q31 HS1 QBl H3 HHI US] is 

an important physical phenomenon occurring inevitably in most experiments 
dealing with quantum objects. It is usually defined as a process whereby the 
physical system of interest interacts with environment or other larger system 
with complex structure and, because of this interaction, changes its evolution 
from unperturbed, coherent internal dynamics. In some sense, the information 
about the initial and subsequent states of system undergoing decoherence is 
leaking into the outer world: The system is no longer described by a wave 
function, but rather by the statistical density matrix [2Ql EO E21 E2 Ei| ■ The 
quantum wave function description only applies to the total system, including 
the environmental modes, which has much more degrees of freedom. Because of 
the importance of quantum coherence for quantum information processing |25[ 

Ha Ha EH Eol EH E21 E31 EH ESI ESI Ell EH ESI Eol EH E21 EH EH EH [66], 

quantitative characterization of decoherence has become an active research 
field with many open problems. 

Since quantum information processing requires maintaining high level of 
coherence, emphasis has recently shifted from large-time system dynamics at 
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experimentally better studied coherence-decay time scales to almost perfectly 
coherent dynamics at much shorter times. Many quantum systems proposed 
as candidates for qubits (quantum bits) for practical realizations of quantum 
computing require quantitative evaluation of their coherence. In other words, a 
single measure characterizing decoherence is desirable for comparison of differ- 
ent qubit designs and their optimization. Besides the evaluation of single qubit 
performance one also has to analyze scaling of decoherence as the register size 
(the number of qubits involved) increases. Direct quantitative calculations of 
decoherence of even few-qubit quantum registers are not feasible. Therefore, 
a practical approach has been to explore quantitative single-parameter mea- 
sures of decoherence [67], develop techniques to calculate such measures at 
least approximately for realistic one- and two-qubit systems [63 [69] , and then 
establish scaling (additivity [70l[7T]) for multi-qubit quantum systems. 

In Section [21 we outline different approaches to define and quantify de- 
coherence. We argue that a measure based on a properly defined norm of 
deviation of the density matrix is appropriate for quantifying decoherence in 
quantum registers. For a semiconductor double quantum dot qubit, evaluation 
of this measure is reviewed in Section [3] For a general class of decoherence 
processes, including those occurring in semiconductor qubits considered in 
Section [U we argue, in Section [U that this measure is additive. Thus, the 
level of quantum noise scales linearly with the number of qubits. 

2 Measures of Decoherence 

In this section, we consider briefly several approaches to quantifying the degree 
of decoherence due to interactions with environment. In Subsection 12.11 we 
discuss the approach based on the asymptotic relaxation time scales. The 
entropy and idempotency-defect measures are reviewed in Subsection l2.2l The 
fidelity measure of decoherence is considered in Subsection 12.31 In Subsection 
12.41 we review our results on the operator norm measures of decoherence. 
Subsection 12 . 51 discusses an approach to eliminate the initial-state dependence 
of the decoherence measures. 

2.1 Relaxation Time Scales 

Decoherence of quantum systems is frequently characterized by the asymp- 
totic rates at which they reach thermal equilibrium at temperature T . One 
of the reasons for focusing on relaxation rates is that large-time behavior is 
relatively easy to observe in ensemble experiments. Markovian approxima- 
tion schemes typically yield exponential approach to the limiting values of 
the density matrix elements for large times [2TJ [22J, [23] . For a two-state sys- 
tem, this defines the time scales T\ and T 2 , associated, respectively, with the 
approach by the diagonal (thermalization) and off-diagonal (dephasing, deco- 
herence) density-matrix elements to their limiting values. More generally, for 
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large times we approximate deviations from stationary values of diagonal and 
off-diagonal density matrix elements as 

Pkk(t) - /Cfcfc(oo) oc e~ t/Tkk , (1) 

p jk (t) ex e"*/ T ^ (j ^ k )- (2) 

The shortest time among Tkk is often identified as T±. Similarly, T2 can be 
defined as the shortest time among T„^ m . These definitions yield the charac- 
teristic times of thermalization and decoherence (dephasing). 

For systems candidate for quantum computing realizations, noise effects 
are commonly reduced by working at very low temperatures and making their 
structure features nanosize for strong quantization. Then for the decoherence 
and thermalization times we have, T2 -C T\, e.g., [21] . Therefore, the decoher- 
ence time is a more crucial parameter for quantum computing considerations. 
The time scale T2 is compared to the "clock" times of quantum control, i.e., 
the quantum gate functions, T g , in order to ensure the fault-tolerant error 
correction criterion T g /T 2 < O (l0~ 4 ), e.g., [66] . 

The disadvantages of this type of analysis are that the exponential be- 
havior of the density matrix elements in the energy basis is applicable only 
for large times, whereas for quantum computing applications, the short-time 
behavior is usually relevant [18]. Moreover, while the energy basis is natural 
for large times, the choice of the preferred basis is not obvious for short and 
intermediate times [HI [72]. Therefore, the time scales T\ and T2 have limited 
applicability in evaluating quantum computing scalability. 



2.2 Quantum Entropy 

An alternative approach is to calculate the entropy [5D] of the system, 

S(t) = -Tr(p\np), (3) 

or the idempotency defect, also termed the first order entropy [73 ] 174 ] [75] . 

S W = l-Tr(p 2 ). (4) 

Both expressions are basis independent, have a minimum at pure states and 
effectively describe the degree of the state's "purity." Any deviation from a 
pure state leads to the deviation from the minimal values, 0, for both measures, 

S pure state(^) — $ pure statc(^) 0- (5) 



2.3 Fidelity 



Writing the total Hamiltonian as follows, 
H = Hs + Hb + Hi , 



(0) 
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where Hg is the term describing internal system dynamics, Hg governs the 
evolution of environment, and Hj describes system-environment interaction, 
let us now define the fidelity [751177]. 

F(t) = Tc s [pui»i{t)p(t)]. (7) 

Here the trace is over the system degrees of freedom, and pidcai(t) represents 
the pure-state evolution of the system under Hs only, without interaction 
with the environment (Hi = 0). In general, the Hamiltonian term Hs govern- 
ing the system dynamics can be time dependent. For the sake of simplicity 
throughout this review we restrict our analysis by constant Hs since approxi- 
mate evaluation of decoherence can be done for qubits controlled by constant 
Hamiltonian. In this case 

Pide a i(t) = e- iHst p(0)e iHst . (8) 

More sophisticated scenarios with qubits evolving under time dependent Hs 
were considered in [78l [79l [80] . 

The fidelity provides a certain measure of decoherence in terms of the dif- 
ference between the "real," environmentally influenced, p(t), evolution and 
the "free" evolution, Pideal(i)- It will attain its maximal value, 1, only pro- 
vided p(t) = Pideal(i)- This property relies on the fact the Pideai(i) remains a 
projection operator (pure state) for all times t > 0. 

As an illustrative example consider a two-level system decaying from the 
excited to ground state, when there is no internal system dynamics, 

/W«)= (9) 



1 



-rt 



P(t) = e -rt ) , (10) 

and the fidelity is a monotonic function of time, 

F(t) = e- n . (11) 

Note that the requirement that pideai(t) is a pure-state (projection opera- 
tor), excludes, in particular, any T > thermalized state as the initial system 
state. For example, let us consider the application of the fidelity measure for 
the infinite-temperature initial state of our two level system. We have 

p(0)=Pi d Ut)=( 1/ Q 2 ^ , (12) 

which is not a projection operator. The spontaneous-decay density matrix is 
then 

P(t)={ o e-r*/2)- (13) 
The fidelity remains constant 
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F(t) = 1/2, (14) 

and it does not provide any information of the time dependence of the decay 
process. 

2.4 Norm of Deviation 

In this subsection we consider the operator norms [81] that measure the devi- 
ation of the system from the ideal state, to quantify the degree of decoherence 
as proposed in [67] . Such measures do not require the initial density matrix 
to be pure-state. We define the deviation according to 

a(t) = p{t) - p ideal (t). (15) 

We can use, for instance, the eigenvalue norm, 

||cr|k = max | Aj| , (16) 

or the trace norm, 

imi Ti . = EM' ( 17 ) 

i 

etc., where Xi are the eigenvalues of the deviation operator (fl5|). A more 
precise definition of the eigenvalue norm for a linear operator, A, is [5T] 



P||=sup 



(18) 



Since density operators are bounded, their norms, as well the norm of the 
deviation, can be always evaluated. Furthermore, since the density operators 
are Hermitian, this definition obviously reduces to the eigenvalue norm ()16p . 
We also note that ||A|| = implies that A — 0. 

The calculation of these norms is sometimes simplified by the observation 
that a(t) is traceless. Specifically, for two-level systems, we get 



roof + h>ir = oll<4&- (is) 



2 , 2 __ -1 

~ 2 

For our example of the two-level system undergoing spontaneous decay, the 



| A = l-e-« (20) 



2.5 Arbitrary Initial States 

The measures considered in the preceding subsections quantify decoherence 
of a system provided its initial state is given. However, this is not always the 
case. In quantum computing, it is impractical to keep track of all the possible 
initial states for each quantum register, that might be needed for implementing 



6 



Leonid Fedichkin and Vladimir Privman 



a particular quantum algorithm. Furthermore, even the preparation of the 
initial state can introduce additional noise. Therefore, for evaluation of fault- 
tolerance (scalability) , it will be necessary to obtain an upper-bound estimate 
of decoherence for an arbitrary initial state. 

To characterize decoherence for an arbitrary initial state, pure or mixed, 
we proposed [67j to use the maximal norm, _D, which is determined as an 
operator norm maximized over all initial density matrices. It is defined as the 
worst case scenario error, 



For realistic two-level systems coupled to various types of environmen- 
tal modes, the expressions of the maximal norm are surprisingly elegant and 
compact. They are usually monotonic and contain no oscillations due to the 
internal system dynamics, as, for example, are the results obtained for semi- 
conductor quantum dot qubits considered in the next section. 

In summary, we have considered several approaches to quantifying deco- 
herence: relaxation times, entropy and fidelity measures, and norms of devia- 
tion, and we defined the maximal measure that is not dependent on the initial 
state, and which will be later shown to be additive; see Section [4j 

3 Decoherence of Double Quantum Dot Charge Qubits 

As a representative example, let us review evaluation of decoherence for semi- 
conductor quantum dots. Quantum devices based on solid-state nanostruc- 
tures have been among the major candidates for large-scale quantum compu- 
tation because they can draw on existing advances in nanotechnology and ma- 
terials processing [82] • Several designs of semiconductor quantum bits (qubits) 
were proposed [23 12H IH3 SH S3 Sll S31 HS1 IHS HU - In particular, the encod- 
ing of quantum information into spatial degrees of freedom of electron placed 
in a quantum dot was considered in [1TJ [321 [531 SS US] ■ A relatively fast decay 
of coherence of electron states in ordinary quantum dots, e.g., [26], can be 
partially suppressed by encoding quantum information in a subspace of elec- 
tron states in specially designed arrays of quantum dots (artificial crystals), 
proposed in [86j . Actually, under certain conditions even double-dot systems 
in semiconductors can be relatively well protected against decoherence due to 
their interactions with phonons and electromagnetic fields [38] . This observa- 
tion was confirmed in recent experiments [45], which demonstrated coherent 
quantum oscillations of an electron in a double-dot structure. 

Several designs of double-dot qubits have been explored in recent experi- 
ments [46j [47l SSI EH [50] carried out at temperatures ranging from tens and 
hundreds of mK. Temperature dependence of relaxation rates in Si charge 
qubits was studied theoretically in [SH [52]. Recently, it has been pointed 
out [68] that in the zero-temperature limit and for conventional double-dot 




(21) 
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structures higher order processes in electron-phonon interaction dominate de- 
coherence. 

In this work, we studied the acoustic phonon bath as the main source of 
decoherence for the considered type of qubit, which is supported by theoretical 
and experimental evidence, e.g., [38l[46]. Decoherence due to different sources, 
e.g., due to trapping center defects [87, 88 , can play important role in other 
situations. 

In the next subsection, we outline the structure of double-dot qubits. Sub- 
sections 13.21 and 13.31 are devoted to the consideration of the electron-phonon 
interaction for two realistic cases: In Subsection l3.2l we analyze the piezoacous- 
tic interaction in crystals with zinc-blende lattice and with parabolic quantum 
dot confinement potential. Double-dots with prevalence of piezointeraction 
have been fabricated [45] in gated GaAs/AlGaAs heterostructures. In Sub- 
section 13.31 we study the deformation interaction with acoustic phonons in 
"quantum dots" formed by double-impurities in semiconductors with inver- 
sion symmetry of elementary lattice cell. Experiments with the latter type of 
double-dot systems have been reported in [471 EH] ■ Finally, Subsections 13.41 
13.51 13.61 and 13.71 present illustrative calculations of the noise level for selected 
quantum gates. 

3.1 Model 

We consider a double-dot structure sketched in Figure [1] It consists of two 
quantum dots coupled to each other via a tunneling barrier and containing 
a single electron hopping between the dots. We limit our consideration to 
double-dot structures in which the energy required to transfer to the upper 
levels is much higher than the lattice temperature and energy spacing between 
the two lowest levels. 



The electron is considered to be in a superposition of two basis states, |0) 
and |1), 



The states that define the "logical" basis are not the physical ground and 
first excited state of the double-dot system. Instead, i/jq (the "0" state of the 
qubit) is chosen to be localized at the first quantum dot and, to a zeroth order 
approximation, be similar to the ground state of that dot if it were isolated. 
Similarly ipi (the "1" state) resembles the ground state of the second dot (if it 
were isolated) . This assumes that the dots are sufficiently (but not necessarily 
exactly) symmetric. We denote the coordinates of the potential minima of the 
dots (dot centers) as vectors Ro and Ri, respectively. The distance between 
the dot centers is 



The Hamiltonian of an electron within a phonon environment is given by 



ip = aip + (3ipi . 



(22) 



L= |L| = |Ri-R |. 



(23) 



H — H e + Hp + H ( 



ep • 



(24) 
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Fig. 1. Electron in a double well potential. 



The electron term is 

H e = ~e A (t)(T x - ^£p(t)a z , (25) 

where a x and a z are Pauli matrices, whereas £^(£) and £p(t) can have time- 
dependent, as determined by unitary single-qubit quantum gate- functions that 
are carried out. They can be controlled externally by adjusting the potential on 
the control electrodes (gates) surrounding the double-dot system. For constant 
ea and ep, the energy splitting between the electron energy levels is 



4^ A 



+ e 2 P . (26) 



The Hamiltonian of the phonon bath is described by 

H p = J2 h ^ b lxKx, (27) 

where x and 6 q t \ are, respectively, the creation and annihilation opera- 
tors of phonons characterized by the wave vector q and polarization A. We 
approximately assume isotropic acoustic phonons, with a linear dispersion, 

uj q = sq, (28) 

where s is the speed of sound in the semiconductor material. In the next 
subsection we show that the electron-phonon interaction can be derived in 
the form 
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with the coupling constants g^ x determined by the architecture of the double- 
dot and the properties of the material crystal structure. 



3.2 Piezoelectric Interaction 

The derivation in this subsection follows 68, 69 . The piezoacoustic electron- 
phonon interaction [89j is described by 

/ * \ 1/2 

H ^= i J2[2^qVj M A(q)^(q)(fo q + ^ q ). (30) 

Here p is the density of the semiconductor material, V is the volume of semi- 
conductor, and for the matrix element M\(q), one can derive 

M A (q) - 2^Ete?j + ^i)%M ijk . (31) 

ijk 

In this expression, ^- are the polarization vector components for polarization 
A, while Mijk express the electric field as a linear response to the stress, 

E k = J2 M ijkSij- (32) 

For a crystal with zinc-blende lattice, exemplified by GaAs, the tensor M^ k 
has only those components non-zero for which all three indexes i, j, k are 
different; furthermore, all these components are equal Mij k = M. Thus, we 
have 

M A (q) = ^-(Ci9 2 93 + ?2<7i93 + ^lfe)- ( 33 ) 

The form factor F(q) accounting for that the electrons in the quantum 
dot geometry are not plane waves, is 

F(q) = Yl 4 C * / d?r^{v)4> k {v)e-^\ (34) 
j,k J 

where c k , Cj are annihilation and creation operators of the basis states k, j — 
0, 1. For gate-engineered quantum dots, we consider the ground states in each 
dot to have an approximately Gaussian shape 

U*) = ^]km^-^^\ (35) 

where 2a is a characteristic size of the dots. 

We assume that the distance between the dots, L, is sufficiently large 
compared to a, to ensure that the different dots wave functions do not overlap 
significantly 



d 3 r0*(r)0 fe (r)e" 



< 1, for j ^ k. (36) 
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This implies that the coupling leading to tunneling between the dots is small, 
as is the case for the recently studied experimental structures [351 SHI S3 EE] , 
where the splitting due to tunneling, measured by ea, was below 20 /xeV, while 
the electron quantization energy in each dot was at least several meV. 
For j = k, we obtain 

2 / 2 



e -*H-n.Jg-« H /■*_ (37) 

The resulting form factor is 

F(q) = e -« 2 « 2 /4 e ^q-R( c t Coet q-L/2 + c t Cie - 4 q-L/ 2)j (38) 

where R = (Ro + Ri) /2. Finally, we get 

= e - a2 9 2 /4 e - 4q -R [ cos ( q . L /2)/ + i S i n (q . L/2)<jJ , (39) 

where / is the identity operator. Only the second term in (|39|) . which is not 
proportional to J, represents an interaction affecting the qubit states. It leads 
to a Hamiltonian term of the form (|29p . with coupling constants 



1/2 



9 ^ = -\^sv) ^- aW4 " iqR 

X (£l e 2 e 3 + ?2 e l e 3 + C3 e l e 2 ) sin (q ' L / 2 ), ( 40 ) 

where e k = q k /q. 

3.3 Deformation Interaction 

Deformation coupling with acoustic phonons [89j is described by 

1/2 



( J^V) ^(^A + ^A), (41) 
q, A 

where S" is a material-dependent constant termed the "deformation potential." 

Here we consider a particular double-dot-like nanostructure which has 
been a focus of recent experiments, due to advances in its fabrication [471 
148) by controlled single-ion implantation: A double-impurity Si structure 
with Hydrogen-like electron confinement potentials for at both impurities (P 
atoms). We consider a Hydrogen- like impurity state, 

^M = ^I72 eHr " R<l/a ' ( 42 ) 

where a is the effective Bohr radius. The form factor in this case is given by 
the following formula, 

F(q) = [l + 6 ( Q V)/4] 2 [C ° S(q ' L/2)/ + * Shl(q ' L/2)CTz] ' (43) 
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The interaction can then be expressed in the form (|29|) . but with different 
coupling constants, 



We note that the details of the electronic band structure, found important 
in recent decoherence studies |90j in other contexts, would matter for calcu- 
lations that either consider decoherence due to processes mediated fully or 
partially by interactions of the bound electron with the conduction electrons, 
or if the impurity-bound electron were hybridized with the conduction elec- 
trons, namely if its wave function were affected by the lattice potential outside 
the dots. In our case, both effects are not considered. Specifically, we assume 
that the quantum dots are strongly binding and that the electron wave func- 
tions are approximately Gaussian or Hydrogen-like. While we do not expect 
the results to be modified significantly, consideration of the band structure 
effects could be an interesting future project especially for P-impurities in 
Si, where the bound-electron wave function actually extends well beyond the 
lattice spacing. 

3.4 Error Estimates During Gate Functions 

In general, the ideal qubit evolution governed by the Hamiltonian term (125|) 
is time dependent. Decoherence estimates for some solid-state systems with 
certain shapes of time dependence of the system Hamiltonian were reported 
recently [751 [79, 80J . However, such calculations are rather complicated. Actu- 
ally, there is no need to consider all possible time dependent controls of qubit 
to evaluate its performance. All single-qubit rotations which are required for 
quantum algorithms can be successfully implemented by using two constant- 
Hamiltonian gates, e.g., amplitude rotation and phase shift [64]. To perform 
both of these gates one can keep the Hamiltonian term ([25]) constant dur- 
ing the implementation of each gate, adjusting the parameters £a and ep as 
appropriate for each gate and for the idling qubit in between gate functions. 

In the following subsections we give specific examples: In Subsection 13.51 
we will consider decoherence during the implementation of the NOT gate (an 
amplitude gate). A 7r-phase shift gate is considered in Subsection 13.61 Then, 
in Subsection 13 . 71 we discuss the overall noise level estimate for a qubit subject 
to gate control. 

3.5 Relaxation During the NOT Gate 

The quantum NOT gate is a unitary operator which transforms the states |0) 
and |1) into each other. Any superposition of |0) and |1) transforms accord- 
ingly, 




1/2 



(44) 



NOT (x|0> + = j/|0> + x|l>. 



(45) 



12 



Leonid Fedichkin and Vladimir Privman 



The NOT gate can be implemented by properly choosing sa and ep in the 
Hamiltonian term (|2l))) . Specifically, with constant 



e A = e (46) 
and 

e P = 0, (47) 

the "ideal" NOT gate function is carried out, with these interaction parame- 
ters, over the time interval T g = r, 

r = — . (48) 

e 

The dominant source of quantum noise for double-dot qubit subject to 
the NOT-gate type coupling, is relaxation involving energy exchange with the 
phonon bath (i.e., emission and absorption of phonons). In this case it is more 
convenient to study the evolution of the density matrix in the energy basis, 
{|+>,|->}, where 



|±) = (|0)±|1))/V2. (49) 

Then, assuming that the time interval of interest is [0, r], the qubit density 
matrix can be expressed |H] as follows, 

fpf + + [ P++ (0) - p%] e~ rt P+ _(0)e-CV 2 -/^ \ 
p(t) = . (50) 

V p- + (0)e-^ r / 2 + is / h ^ f&_ + [p— (0) - p^_] e- rt J 

This is the standard Markovian approximation for the evolution of the density 
matrix. For large times, this evolution would result in the thermal state, with 
the off-diagonal density matrix elements decaying to zero, while the diagonal 
ones approaching the thermal values proportional to the Boltzmann factors 
corresponding to the energies ±e/2. However, we are only interested in such 
evolution for a short time interval, r, of a NOT gate. The rate parameter r 
is the sum [22] of the phonon emission rate, W e , and absorption rate, W a , 

r = w e + w a . (51) 

The probability for the absorption of a phonon due to excitation from the 
ground state to the upper level is 
2tt, 

where \i) is the initial state with the extra phonon with energy fisq and |/) is 
the final state, q is the wave vector, and A is the phonon polarization. Thus, 
we have to calculate 



w X = — \(f\H ep \i)\ 2 5(e-hsq), (52) 



W a = > w 



^ (2tt) 3 



£ / d 3 qw x . (53) 



For the interaction ([29)) one can derive 
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2tt, 



10 : 

where 
N th 



1 



(54) 



(55) 



exp(hsq/kgT) — 1 

is the phonon occupation number at temperature T, and Ub is the Boltzmann 
constant. 

For the piezoacoustic interaction, the coupling constant in (|40p depends 
on the polarization. For longitudinal phonons, the polarization vector has 
Cartesian components, expressed in terms of the spherical-coordinate angles, 



£i =e.\ = sin 9 cos 4>, (£ — e 2 = sin 9 sin 0, £ 3 — e 3 = cos0, 



(56) 



where ej = qj/q. For transverse phonons, it is convenient to define the two 
polarization vectors £j 1 and ^ 2 to have 



Cf 1 - sin 
2 



COS COS ( 



c ±2 
S2 



cos f sin ( 



sm( 



Then for longitudinal phonons, one obtains [69] 

w ll = — — — Af 2 e — a2 * 2 ^ 4 
psVq 

x 9 sin 4 9 cos 2 sin 2 </> cos 2 (j) sin 2 (gL cos 9/2) . 
For transverse phonons, one gets 



psVq 



M 2 e -aV/4(_ 2 



sin cos W sin cos < 



w 



sin 3 #cos0sin <ft) 2 sin 2 (gL cos#/2), 
77 M 2 e~ a2<z2/4 (-2sin6lcos6'cos 2 ( / 



+ sin cos 9 sin 2 <^>) 2 sin 2 (gi cos 9/2) . 

By combining these contributions and substituting them in 
probability of absorption of a phonon for all polarizations, 



M 2 ex P 
W a - = — 

pl0ZO 20-Kps 2 hL 5 k i 



exp 



{(kL) 1 



( risk \ 

\k B T ) 



- 1 



5kL 



2 (kL) - 21 cos (kL) 
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3(kL) 2 sin(kL)X, 



where 



(57) 
(58) 

(59) 



(60) 



(61) 
we get the 



(62) 
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is the wave-vector of the absorbed phonon. 

For the deformation interaction one can obtain the following result, 

W = ^ [l + (a V)/4] 4 Sin2 ( q ' ~ ™ 

The total probability for a phonon absorption is 

_ H 2 fc 3 l-sin(fc£)/(fcL) 

"deform ~ / i„„„2fc 71 o, n ,,»4 / I , \ ' V D0 7 



1 



Finally, the expressions for the phonon emission rates, W e , can be obtained 
by multiplying the above expressions, (|62j) and (f65|) . by (JV t /, + 1)/N t h- 



3.6 Dephasing During a Phase Gate 

The 7r gate is a unitary operator which does not change the absolute values of 
the probability amplitudes of a qubit in the superposition of the |0) and |1) 
basis states. It changes the relative phase between the probability amplitudes. 
Specifically, any superposition of |0) and |1) transforms according to 

/T(x|0)+y|l»=a;|0>-y|l>. (66) 

Over a time interval r, the 7r gate can be carried out with constant interaction 
parameters, 

£A = (67) 
and 

irh 

e P =e = — . (68) 

T 

In [55] , double-dot qubit dynamics during implementation of phase gates 
was considered. The relaxation dynamics is suppressed during the ir gate, 
because there is no tunneling between the dots. Quantum noise then results 
due to pure dephasing, i.e., via the decay of the off-diagonal qubit density 
matrix elements, while the diagonal density matrix elements remain constant. 
In the regime of pure dephasing, the qubit density matrix can be represented 

as mm] 

/ Poo(0) P01 (0)e- B2(t)+tEt/h \ 

P (t) = , (69) 

\p 1Q (0)e- B2 ^-^/ h Pll (0) / 

with the spectral function, 

B \t) = \y W! sin 2 coth 4^ 

W h 2 ^ uj 2 „ 2 2k B T 

q,A y 

" ' (70) 



»V J q's' 2 2k B T' 
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For the piezoelectric interaction, the coupling constant \ was obtained 
in (jit)]) , and expression for the spectral function takes the form 

POO P7T /"27T 



x ^ (^a+^+^i^) 2 exp(-aV /2) 

A q 



x sin (gL cos ft) sin — coth — , (71) 

2 IkbT 

c.f. ([56)) - ([58]) . For the deformation interaction, we have the coupling constant 
(|44[) , and the expression for the spectral function is given by 

=•2 



S deform(*) = 2 T 3 / 9 2rf 9 / sin 9de 
7T Alps Jo JO 



sin 2 (oi cos #) . o osi S,gs 



3.7 Qubit Error Estimates 



The qubit error measure, D, is obtained from the density matrix deviation 
from the "ideal" evolution by using the operator norm approach |67j reviewed 
in Subsection 12.51 After lengthy intermediate calculations one gets j68] rela- 
tively simple expressions for the error during the NOT gate, 

1 - e- rT 

^NOT = ; , _ e/kBT , (73) 



1 + e 
and the 7r gate, 
1 
2 



A, 



l- e -B 2 (r) 



(74) 



Table 1. Qubit parameters 



Parameter 


GaAs double-dot qubit 


Si double-impurity qubit 


p, kg/m 3 


5.31 x 10 3 


2.33 x 10 3 


s, m/s 


5.14 x 10 3 


9.0 x 10 3 


S, eV 


3.3 




ei4, C/m 2 




0.16 


«: 




12.8 


M, eV/m 




eei4/(eo«) 


L, nm 


50 


50 


a, nm 


25 


3 
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Fig. 2. Estimates of the error measure per cycle, D, due to the piezoelectric inter- 
action in GaAs double-dot, shown as a function of the temperature, T. The cycle 
time t was 6 • 10 -11 s. 



A realistic noise estimate could be taken as the worst case scenario, i.e., the 
maximum of these two expressions for error per gate cycle. The expressions 
(|75|) and (f?4")) were used to calculate the error rate for the double-dot qubit in 
GaAs and double-impurity qubit in Si. The parameters used were chosen to 
correspond to the experimentally realized structures, [46j |45l [47l [48] , and are 
summarized in Table 1 . The calculated error measures are presented in Figures 
2 and 3. The gate time r selected for the reported calculations, 6- 10 -11 s, is a 
representative value consistent with typical experimental conditions. In fact, 
decreasing the gating time does not lead to smaller quantum noise in this case 
because the energy gap of the driven qubit is ~ 1/r. If the gap is made too 
large, other excitations will play a role in decoherence, for instance, optical 
phonons. The time scale chosen here is within an optimal range, as discussed 
in [55]. 

In summary, we derived expressions for the error measure for double-dot 
and double-impurity qubits. The results, presented in Figures 2 and 3, sug- 
gest that pure dephasing dominates at low temperatures. As the temperature 
increases beyond about 1 K, the effect of relaxation becomes comparable and 
ultimately dominant. 
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Fig. 3. Estimate of error rate per cycle, D, due to deformation phonon interaction 
for a double Phosphorus impurity in Si, shown as a function of the temperature, T. 
The cycle time r was 6 • 10~ n s. The relaxation rate for this range of the parameter 
values is negligibly small and respective values of D are not shown. 

The error measure values found, are 1.5 or more orders of magnitude larger 
than the "traditional" fault-tolerance thresholds for multiqubit quantum com- 
putation, which range from O(10- 4 ) down to O(10- 6 ) [5J M EH Mi E>1 E3] • 
However, recent developments have yielded less strict requirements for the 
error rate [Ml [Ml HIS], optimistically, as large as O(10 -2 ). Furthermore, 
there are several approaches to decrease decoherence effects by pulsed control 
[971 [23 [Ml H001 ESI EU EH EH1 Ell EU, some recently tested exper- 
imentally in multi-spin NMR [107, 108 . Other ideas rely on the fact that 
instead of the bulk material, the qubit could be manufactured in a one- or 
two-dimensional nanostructure [1091 1110"] . the latter already available experi- 
mentally jl 1 1] , which would affect the phonon spectrum and lower decoherence 
effects. 



4 Additivity of Decoherence Measures 

In the study of decoherence of several-qubit systems, additional physical ef- 
fects should be taken into account. Specifically, one has to consider the degree 
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to which noisy environments of different qubits are correlated [911 1112) . In ad- 
dition to acting as a source of the quantum noise, the correlated bath can in- 
duce an effective interaction, namely, create entanglement, between the qubits 
immersed in it [1101 11131 1114L ITT5] . Furthermore, if all constituent qubits are 
effectively immersed in the same bath, then there are ways to reduce decoher- 
ence for this group of qubits without error correction algorithms, by encoding 
the state of one logical qubit in a decoherence-free subspace of the states of 
several physical qubits [86l [HI 11161 1117) 1118] . In this section, we will consider 
several-qubit quantum registers and, as the "worst case scenario" assume that 
the qubits experience uncorrelated noise, i.e., each is coupled to a separate 
bath. Since analytical calculations for several qubits are not feasible, we seek 
"additivity" properties that will allow us to estimate the error measure for 
the register from the error measures of the constituent qubits. 

It is important to emphasize that loss of quantum coherence results in a 
loss of various two- and several-qubit entanglements in the system. The high- 
est order (multi-qubit) entanglements are "encoded" in the far off-diagonal 
elements of the multi-qubit register density matrix, and therefore these quan- 
tum correlations will decay at least as fast as the products of the decay fac- 
tors for the qubits involved, as exemplified by several explicit calculations 
|1191 1120"! 11211 [T2"2"] . This observation leads to the conclusion that, for large 
times, the rates of decay of coherence of the qubits will be additive. 

However, here we seek a different result: one valid not in the regime of the 
asymptotic large-time decay of quantum coherence, but for relatively short 
times, t, of quantum gate functions, when the noise level, namely the value 
of the measure D(r) for each qubit, is relatively small. In this regime, we will 
establish [70j in this section, that, even for strongly entangled qubits — which 
is important for the utilization of the power of quantum computation — the 
error measures D of the individual qubits in a quantum register are additive. 
Thus, the error measure for a register made of similar qubits, scales up lin- 
early with their number, consistent with other theoretical and experimental 
observations [H [Ml CEDE] . 

In Subsection 14. 1[ we revisit the noise measure via the maximal deviation 
norm and discuss some of its properties. In Subsection 14. 2\ we introduce the 
diamond norm which is used as an auxiliary tool in the proof of additivity. We 
then establish an approximate upper bound for D(t) for a register of several 
weakly interacting but possibly strongly entangled qubits, and cite work that 
further refines the additivity properties for typical qubit realizations. 

4.1 The Maximal Deviation Norm 

To characterize decoherence for an arbitrary initial state, pure or mixed, we 
use the maximal norm, D, which was defined (|2 1 j) in Subsection 12.51 as an 
operator norm maximized over all the possible initial density matrices. One 
can show that < D(i) < 1. This measure of decoherence will typically in- 
crease monotonically from zero at t — 0, saturating at large times at a value 
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D (oo) < 1. The definition of the maximal decoherence measure D(t) looks 
rather complicated for a general multiqubit system. However, it can be evalu- 
ated in closed form for short times, appropriate for quantum computing, for a 
single-qubit (two-state) system. We then establish an approximate additivity 
that allows us to estimate D(t) for several-qubit systems as well. 

In the superoperator notation the evolution of the reduced density operator 
of the system ([7]) and the one for the ideal density matrix ((8]) can be formally 
expressed [62l |63l [64] in the following way 

p(t) = T(i)p(Q), (75) 



pW( t ) =T W(t)p(0), 



(76) 



where T, T« are linear superoperators. In this notation the deviation can be 
expressed as 

a(t)=[T(t)-T^(t)]p(0). (77) 
The initial density matrix can always be written in the following form, 



P(0) 



E 



(78) 



where J2jPj = 1 an d < pj < 1. Here the set of the wavefunctions \ipj) is 
not assumed to have any orthogonality properties. Then, we get 



a(i, P (o)) = 5>, [r(t)-rW(*)] l^x^-l. 

3 



The deviation norm can thus be bounded. 



W(t,p(o))\\ 



< 



r(t)-rW(t)l |0X0| 



Here |0) is defined according to 



rp _ Jl(j) 



max 



(79) 



(80) 



(81) 



It transpires that for any initial density operator which is a statistical mixture, 
one can always find a density operator which is pure-state, \<fi){(t>\, such that 
\\cr(t, p(0))\\x < || (r(t, |0) (0|) || a ■ Therefore, evaluation of the supremum over 
the initial density operators in order to find D(t), see (12TT) , can be done over 
only pure-state density operators, p(0). 

Let us consider strategies of evaluating D(t) for a single qubit. We can 
parameterize p(0) as 



p(0) = U 



P 
1-P 



where < P < 1, and [/ is an arbitrary 2x2 unitary matrix, 
' e i ( Q +7) C0S 6) e ^ a -^ sin 6 



U 



a i(j-a) 



sm( 



-i(a+7) 



COS ( 



(82) 



(83) 
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Then, one should find a supremum of the norm of deviation (|16|) over all the 
possible real parameters P, a, 7 and 8. As shown above, it suffices to consider 
the density operator in the form of a projector and put P = 1. Thus, one 
should search for the maximum over the remaining three real parameters a, 
7 and 9. 

Another parametrization of the pure-state density operators, p(0) = 
\(t>)(4>\ : is to express an arbitrary wave function \<fi) = J2j( a j m some 

convenient orthonormal basis where j = 1, . . . , N. For a two-level system, 



where the four real parameters 01,2,61,2 satisfy a\ + b\ + a\ + b\ = 1, so 
that the maximization is again over three independent real numbers. The 
final expressions ([73]) and ([74)1 for D(t), for our selected single-qubit systems 
considered in Section [31 are actually quite compact and tractable. 

In quantum computing, the error rates can be significantly reduced by us- 
ing several physical qubits to encode each logical qubit [HS1 11161 1117] . There- 
fore, even before active quantum error correction is incorporated [531 (Ml [551 
l56l l57l l58l l59l IBTJl |6T] , evaluation of decoherence of several qubits is an impor- 
tant, but formidable task. Thus, our aim is to prove the approximate additivity 
of D q {t), including the case of the initially strongly entangled qubits, labeled 
by whose dynamics is governed by 



where Hg q is the Hamiltonian of the qth qubit itself, H^q is the Hamiltonian of 
the environment of the qth qubit, and Hj q is corresponding qubit-environmcnt 
interaction. In the next subsection we consider a more complicated (for actual 
evaluation) diamond norm 62, 63, 64 , K(t), as an auxiliary quantity used to 
establish the additivity of the more easily calculable operator norm D(t). 

4.2 Upper Bound for Measure of Decoherence 

The establishment of the upper-bound estimate for the maximal deviation 
norm of a multiqubit system, involves several steps. We derive a bound for 
this norm in terms of the recently introduced (in the context of quantum 
computing) 62, 63, 64 diamond norm, K(t). Actually, for single qubits, in 
several models the diamond norm can be expressed via the corresponding 
maximal deviation norm. At the same time, the diamond norm for the whole 
quantum system is bounded by sum of the norms of the constituent qubits 
by using a specific stability property of the diamond norm. The use of the 
diamond norm was proposed in [62, 63, 64 , 




(84) 




(85) 



K(t) = \\T- T<% = sup || {[T - T^}®I}g\\ Tr . 



(86) 



Q 
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The superoperators T, TW characterize the actual and ideal evolutions ac- 
cording to [[75]). (fTS 7 )) . Here / is the identity superoperator in a Hilbert space 
G whose dimension is the same as that of the corresponding space of the su- 
peroperators T and T^ % > , and g is an arbitrary density operator in the product 
space of twice the number of qubits. 

The diamond norm has an important stability property, proved in [62[ I63|, 

EH- 

||Bi*B 2 |U = ||BiW|Ba||o. ( 8? ) 



Note that (|87|) is a property of the superoperators rather than that of the 
operators. 

Consider a composite system consisting of the two subsystems Si, S2, with 
the noninteracting Hamiltonian 



iSi i" Hs 2 . 

The evolution superoperator of the system will be 

T Sl s 2 =T Sl ®T S2 , (89) 
and the ideal one 

tig* (90) 

The diamond measure for the system can be expressed as 

K Sl s 2 = \\T Sl s 2 -T { s% 2 \\» = \\(T Sl -T ( s!>T S2 +T^^T S2 -T ( i% 

< \\(t Si - t£)®t S2 \u + \\t£ht S2 - (91) 

By using the stability property (|57j). we get 

K SlS2 < \\(T Sl T^T S2 \\ + \\T^(T S2 - T$>)\U - 
\\T Sl -T^UT s J\. + \\T^%\\T S2 -T^\U = 

\\T Sl - 4% + \\T S2 - iglU = K Si + K S2 . (92) 



The inequality 



K<Y,K q , (93) 



for the diamond norm K(t) has thus been obtained. Let us emphasize that 
the subsystems can be initially entangled. This property is particularly useful 
for quantum computing, the power of which is based on qubit entanglement. 
However, even in the simplest case of the diamond norm of one qubit, the 
calculations are extremely cumbersome. Therefore, the use of the measure 
D(i) is preferable for actual calculations. 

For short times, of quantum gate functions, we can use (|9"3"|) as an approx- 
imate inequality for order of magnitude estimates of decoherence measures, 
even when the qubits are interacting. Indeed, for short times, the interaction 
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effects will not modify the quantities entering both sides significantly. The 
key point is that while the interaction effects are small, this inequality can be 
used for strongly entangled qubits. 

The two deviation-operator norms considered are related by the following 
inequality 

Ha < I \WK < !■ (94) 

Here the left-hand side follows from 

Tra = ^A j= 0. (95) 
j 

Therefore the £th eigenvalue of the deviation operator a that has the maximum 
absolute value, Xi — A ma x, can be expressed as 



Thus, we have 



The right-hand side of (|94|) then also follows, because any density matrix has 
trace norm 1, 

NlTr = \\P - p {l) \Wi < WpWtt + ||p W ||Tr = 2. (98) 

From the relation (|98p it follows that 

K(t) < 2. (99) 
By taking the supremum of both sides of the relation (|97|) we get 

D(t) = sup |H| A < \ sup ||a|| Tr < \K{t), (100) 
p(o) z p(0) z 

where the last step involves technical derivation details [7D] not reproduced 
here. In fact, for a single qubit, calculations for typical models [70] give 

1 

Since D is generally bounded by (or equal to) K/2, it follows that the mul- 
tiqubit norm D is approximately bounded from above by the sum of the 
single-qubit norms even for the initially entangled qubits, 



D q {t) = Mt). (101) 



D(t) < l -K{t) < \ ^ K q (t) = Y, DM (102) 
where q labels the qubits. 
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For specific models of decoherence of the type encountered in Section [31 as 
well as those formulated for general studies of short-time decoherence [57], a 
stronger property has been demonstrated [70] , namely that the noise measures 
are actually equal, for low levels of noise, 



In summary, in this section we considered the maximal operator norm suit- 
able for evaluation of decoherence for a quantum register consisting of qubits 
immersed in noisy environments. We established the additivity property of 
this measure of decoherence for multi-qubit registers at short times, for which 
the level of quantum noise is low, and the qubit-qubit interaction effects are 
small, but without any limitation on the initial entanglement of the qubit 
register. 
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